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ABSTRACT 


Creep  data  of  2618-T61  aluminum  alloy  under  multistep  multiaxial 
proportional  loadings  at  200°C  (392°F)  are  reported.  Two  viscoplastic 
flow  rules  were  developed  using  constant  stress  creep  and  strain 
recovery  data.  One  was  based  on  the  accumulated  strain  (strain 
hardening),  and  the  other  on  a  tensorial  state  variable  (kinematic 
hardening).  Data  were  represented  by  two  models:  a  nonrecoverable 
viscoplastic  model;  and  a  viscous-viscoelastic  model  in  which  the 
time-dependent  strain  was  resolved  inLo  recoverable  (viscoelastic)  and 
nonrecoverable  components.  The  modified  superposition  principle  was 
used  to  predict  the  viscoelastic  strain  component  under  variable 
stress  states  for  both  models.  The  experiments  showed  that  the 
viscous-viscoelastic  model  with  either  strain-hardening  or  kinematic 
hardening  gave  very  good  predictions  of  the  material  responses. 

Strain  hardening  was  best  in  some  step-down  stress  states.  The 
viscoelastic  component  accounted  for  not  only  the  recovery  strain  but 
also  the  transient  creep  strain  upon  reloadings  and  step-up  loadings. 


INTRODUCTION 


There  have  been  great  advances  in  development  of  constitutive 
relations  suitable  for  material  used  in  design  for  high  temperature. 
However,  all  suffer  from  a  lack  of  experimental  data  to  evaluate  them 

[1J. 

Classically,  most  constitutive  relations  were  based  on  the  strain 
hardening  rule.  Some  modifications  have  to  be  made  for  this  approach 
to  simulate  the  anisotropic  nature  observed  in  most  structural 
materials,  see  [2,3]  and  [4,  Chap.  2]. 

In  recent  theories  state  variables  have  been  introduced  to 
incorporate  the  anisotropic  nature.  The  formulations  usually  consist 
of  a  flow  rule  which  describes  the  interrelations  between  strain  rate, 
applied  stress,  and  state  variables,  and  an  evolution  rule  which 
describes  the  rate  of  change  of  the  state  variables  with  deformation 
and  metallurgical  effects.  The  evolution  rule  was  generally  based  on 
the  Bailey-Orowan  model  [5],  and  in  most  cases  a  tensorial  state 
variable  (back  stress)  was  associated  with  the  kinematic  hardening 
nature  and  a  scalar  state  variable  with  the  isotropic  hardening 
nature,  see  [6-8] .  Nearly  all  theories  used  the  Mises  relation  to 
extend  their  applicability  from  uniaxial  to  multiaxial  stress  states, 
and  except  for  Hart's  work  [8],  they  are  all  based  on  the  viscoplastic 
approach. 

In  previous  work  at  Brown  University  on  2618-T61  aluminum  [9-12], 
a  viscous-viscoelastic  model  (W)  with  strain  hardening  nature  was 
developed  and  compared  with  a  series  of  short  term  creep  tests  under 
various  stress  histories.  The  stress  strain  relations  were  refined  in 
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[13]  through  a  set  of  48  hour  creep  and  recovery  data  under  constant 
multiaxial  stress  states  including  aging  and  compression  tests.  In 
current  work,  the  previous  strain  hardening  flow  rule  (SH)  was  refined 
and  another  kinematic  hardening  rule  (KH)  was  derived  from  the  stress 
strain  relations  reported  in  [l3j.  In  addition  to  the  W  model,  a 
viscoplastic  model  (VP)  in  which  all  the  time  dependent  strain  was 
assumed  nonrecoverable  was  also  developed  to  show  the  role  of  the 

i 

viscoelastic  component.  A  set  of  300-hour  creep  experiments  under 
multiaxial  multistep  proportional  loadings  were  compared  with:  a  W 
model  with  SH  nature  (W-SH);  a  W  model  with  KH  nature  (W-KH);  a  VP 
model  with  SH  nature  (VP-SH);  and  a  VP  model  with  KH  nature  (VP-KH). 
Future  work  will  consider  more  complicated  stress  states  such  as 
nonproportional  loadings,  stress  reversals  and  mixed  stress  relaxation 
and  multiaxial  creep. 

MATERIAL,  SPECIMEN,  APPARATUS  AND  PROCEDURES 
The  material  employed  in  the  present  work  was  aluminum  alloy 
2618-T61  which  came  from  the  same  batch  as  that  reported  in  [ 13] ,  but 
probably  different  from  that  used  in  [9-12].  A  more  complete 
description  of  material  and  specimen  were  given  in  [9,13].  The 
combined  tension  and  torsion  machine  and  the  compression  machine  used 
for  these  experiments  were  described  in  [14]  and  [15],  respectively. 
The  temperature  control  system  and  the  data  acquisition  procedures 
were  described  in  [9,14,15].  All  the  tests  were  performed  at  a  test 
temperature  of  200°C. 
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EXPERIMENTAL  RESULTS 

Eight  tests  are  reported  here  whose  loading  programs  and  the 
resulting  total  strain  versus  time  are  shown  in  Fig.  1-8.  The  total 
duration  of  each  test  was  around  300  h.  The  letter  A  preceeding  the 
test  number  represents  a  pure  tension  or  compression  test,  T 
represents  a  pure  torsion  test,  and  CA,  CT  represent  the  axial  part 
and  the  torsion  part,  respectively,  in  a  combined  tension  and  torsion 
test.  The  first  two  steps  of  each  test  were  constant  stress  creep 
followed  by  strain  recovery.  These  were  used  for  determining  the 
stress-strain-time  relations  for  constant  stress  states,  as  reported 
in  [13].  Test  23,  Fig.  1,  was  not  shown  in  [13]  because  some  friction 
was  detected  in  the  strain  measuring  system  during  the  first  step  and 
small  adjustments  were  made  during  step  2  to  save  the  following  steps. 

The  loading  direction  remains  constant  for  each  test.  The  stress 
magnitudes  change  abruptly  in  the  following  ways:  step  up  followed  by 
step  down  as  in  Fig.  2,  step  down  followed  by  step  up  as  in  Fig.  3  and 
8,  and  step  down  followed  by  reloading  and  step  down  as  in  Fig.  1, 

4-7.  Test  A37,  Fig.  7,  is  a  compression  test  with  the  loading  program 
exactly  the  same  as  that  of  tension  test  A32  shown  in  the  same  figure. 
The  symmetry  in  tension  and  compression  as  discussed  in  [13]  can  be 
reaffirmed  from  Fig.  7.  Another  compression  test  is  A44,  Fig.  8. 

The  loading  program  for  test  31  (CA31  or  CT31),  Fig.  5,  6  is 
basically  the  same  as  that  for  test  A32  with  equal  magnitudes  of 
effective  stress  defined  by  *  3x^  ■  (°eff^  *or  eac**  P*ir  of 
corresponding  steps. 

The  third  step  for  all  the  tests  except  test  A24  Fig.  2  was 
reloaded  to  the  same  stress  as  the  first  step.  This  always  resulted 
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in  some  transient  creep  strain  upon  reloading  to  the  previous  creep 
stress  level.  Also  all  the  step  up  steps  resulted  in  transient  creep. 
Partial  unloading  resulted  in  the  following:  If  the  reduction  was 
small  enough  there  was  a  period  of  strain  recovery  followed  by  new 
creep  within  24  h  or  less,  as  in  period  6  Fig.  1  and  2,  period  4  Fig.  3, 
period  7  Fig.  5,  6  and  7.  However,  in  period  4  Fig.  8  a  reduction  of 
stress  similar  to  that  of  period  4  Fig.  3  but  at  much  higher  stress 
yielded  continuous  strain  recovery  during  the  period  of  observation. 

When  the  stress  was  reduced  in  one  or  two  steps  to  about  one  half 
there  was  a  short  strain  recovery  period  followed  by  a  long  period 
with  no  strain  change,  as  in  period  7  Fig.  1  and  2,  period  4  Fig.  4, 

5,  6,  and  7,  period  7  Fig.  4,  and  period  9  Fig.  8.  If  the  stress  was 
reduced  in  several  steps  to  a  sufficiently  small  proportion  of  the 
initial  stress,  continuous  slow  strain  recovery  was  observed  as  in 
period  8  Fig.  1,  6  and  7,  period  5  Fig.  3  and  8. 


STRESS-STRAIN-TIME  RELATIONS  FOR  CONSTANT  STRESS  STATES 
In  the  W  model,  the  total  strain  e^j  during  constant  stress  creep 
was  represented  as: 


c  •  • 


E  +VE  "l  +V  n2 
eij  ♦  eij  t  +  Hjt  , 


(l) 


E  .  .  .  .  +VE  "l  VE 

where  eij  is  the  time  independent  elastic  strain,  eij  t  ■  eij 

+V  n2  V 

viscoelastic  component,  and  cj,jt  ■  eij  is  the  time  dependent 


is  the 


nonrecoverable  component  for  the  W  model,  in'  we  0.223  and 

0.496,  respectively,  see  [13].  Time  independent  plastic  strain  was 
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negligible  for  this  material  at  the  stresses  employed.  In  the  VP  model, 
the  total  strain  during  creep  is 


•  +  p*  •  t-N 

eij  *  E1JC 


(2) 


+  y 

where  e^jt^  *  e^j  is  the  time  dependent  strain  for  the  VP  model,  and  N 
was  found  to  be  0.407,  [ 13] .  The  elasLic  moduli  for  tension  E  and  for 
torsion  G  were  found  to  be  6.50x10^  MPa  (9.43x10^  ksi)  and  2.38x10^  MPa 
(3.45x10^  ksi),  respectively.  The  third  order  multiple  integral 


representation  yields  the  following  expressions  for  the  stress  dependent 
♦VK  +VE 

coefficients  and  e^2  under  constant  combined  tension  and  torsion 

stress,  i.e., 


F^O,!) 


e12F  “  Gve(o,t)  *  Gji  +  G^t^  +  G4(j2t  ,  (4) 


where  F£  and  Gj  are  constants  whose  values  are  given  in  [13] . 

A  function  Fa  (or  Ffc)of  maximum  shear  stress  xmnv.  multiplied  by 

the  stress  deviator  0£j,  was  developed  to  represent  the  stress 
♦V  + 

dependence  of  c^j  or  C£j  as  follows, 

♦V  ' 

eij  *  ®ijEa^Tmax^  » 

where 

Fa( Tmax)  “  3.616xl0“5x{l+0.3914x  exp[2. 108xl05x(Tmax/G)2] }  , 

percent  per  MPa-hr^*^®®  , 
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4 


I 


where  G  is  the  shear  modulus.  Equations  (5)  and  (6)  represent  the 
data  best  of  several  different  formulations  considered,  such  as  power 
function,  hyperbolic  sine  and  combinations  of  power  and  exponential 
functions. 

WEAKENING  EFFECTS 

In  previous  work  [13],  secondary  and  tertiary  creep  stages  were 
observed  when  the  magnitude  of  the  effective  stress  was  137.2  MPa 
(19.90  ksi)  or  higher.  It  was  also  reported  that  there  were  no 
significant  aging  effects  in  this  precipitation  hardened  material  in 
300  hours.  The  onset  of  the  secondary  and  tertiary  stages  could  be 
caused  by  a  combination  of  several  other  weakening  effects  such  as 
static  recovery,  dynamic  recovery,  grain  boundary  shearing  and 
intercrystalline  fracture.  In  current  work,  these  possible  weakening 
effects  were  neglected  in  the  theoretical  analysis,  which  did  not 
consider  the  occurrence  of  secondary  and  tertiary  stages. 
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CONSTITUTIVE  EQUATIONS  FOR  VARIABLE  STRESS  STATES 
VE 

VISCOELASTIC  STRAIN,  e£j:  By  the  modified  superposition  principle 
[9,  16,  17]  the  viscoelastic  strain  for  a  varying  stress  is  given  by 


e 


I 


8f [o(S).t-g] 
3o(£) 


a(OdC 


(7) 


where  the  strain  at  constant  stress  is  given  by  e  =  f(o,t). 

For  a  series  of  m  steps  in  stress  as  in  the  present  test  programs, 
the  strain  was  described  by  the  following  form  from  Eq.  (7) 


-  FVECoijTxHt  1  -  (t-ti)  1]  +  ... 

nl  nl 

+  FVE(om_i,T,n«i)[(t-tni-2^  “  (t-tjn-i)  ] 

nl 

+  F^WnKt-to-i)  ,  tn-i  <  t  .  (8) 


Similarly,  £\2  was  obtained  by  replacing  F^®(o,x)  by  G^®(o,t)  in 
Eq.  (8). 


V 

TIME  DEPENDENT  NONRECOVERABLE  STRAIN,  £ij:  Two  flow  rules, 

V 

strain  hardening  or  kinematic  hardening,  were  developed  here  for 
under  variable  stress  states.  The  derivation  was  based  on  the 
assumption  that  the  hardening  nature  of  the  material  under  variable 
stress  states  was  the  same  as  that  under  constant  stress  states. 


DEVELOPMENT  OF  THE  STRAIN  HARDENING  FLOW  RULE:  Here  only  the 
V 

case  for  the  ejj  (VP  model)  is  discussed.  The  derivations  are  exactly 
V 

the  same  for  the  (W  model)  except  for  different  values  for  the 

time  exponents  and  different  relationships  for  F(Tmax). 

V 

As  shown  in  [13],  &ij  under  constant  stress  creep  can  be 
expressed  as: 

cij  “  F<Tmax>aijtN  •  (9) 


Differentiating  both  sides  with  respect  to  time  yields 


NF(Tmax 


(10) 


To  obtain  a  strain  hardening  flow  rule  the  t  in  Eq.  (10)  was 
expressed  in  terms  of  the  accumulated  strain  as  follows.  Taking  the 
maximum  of  the  shear  strain  on  the  left  hand  side  of  Eq.  (9)  and  the 
|  maximum  of  the  shear  stress  on  the  right  hand  side  yields  the 

following  expression: 

eij  “  Tmax^(Tmax)t^  » 

I 


from  which, 
t 


t e  i  j  /  [  Tma::F  ( Tmax)  ]  ^ 


1/N 


(12) 


where  Ejj  is  the  magnitude  of  the  strain  tensor  defined  by  the  Tresca 


relation.  The  choice  of  the  Tresca  relation  instead  of  a  Mises 
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relation  for  magnitude  quantities  came  from  the  observation  that  the 
maximian  shear  rates  are  constant  at  any  given  time  for  those  creep 
tests  whose  stress  levels  are  the  same  according  to  the  Tresca  relation 
as  discussed  in  [13] .  For  combined  tension  a  and  torsion  t  stress 
states,  Tmav  is  equal  to  [(a/ 2)^+t2]1/2  and  ejj  equals 
{[(3/4)eii]2+(ei2^]^^*  Substituting  t  from  Eq.  (12)  into  Eq.  (10) 
yielded  the  strain  hardening  flow  rule, 


where  F(Tmj,T)  is  as  follows:  For  the  VP  model,  F(Tmax)  is  F^t,,^), 
Eq.  (6).  For  the  W  model,  N  is  replaced  by  n2  and  F(Tmax)  is 
Fa  (T  max) »  Eq.  (5). 

Furthermore,  differentiating  both  sides  of  Eq.  (11)  with  respect 

V 

to  time  t,  substituting  Eq.  (12)  for  t  and  dividing  by  ejj  yields 


V  V 

eii/t (ei ' 


} (N-l) /Nj 


Integrating  (14)  with  respect  to  time  yields 

-  {h*max<OF(TmaxU))]l/V}N  .  (15) 

0 

For  a  series  of  m  steps  in  stress,  Eq.  (15)  becomes 

~  l/H 

eij  *  ( [Tmax(l)^(Tmax(l))]  <fcl)  +  ••• 

1/N 

+  tTmax(m-l)^(',max(m-l))l  (tm-l“tm-2) 

r  « 1/N  N 

+  lTmax(m)E(Tmax(m))J  (t-t^-i)) 


,  ta-1  <  t  . 


(16) 
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Eq.  (16)  is  true  for  both  proportional  and  nonproportional  loadings. 
However,  for  proportional  loadings,  the  strain  components  at  any  given 
time  can  be  obtained  directly  from  Eq.  (16)  because  the  loading  direction 
is  constant  throughout  the  test.  From  Eq.  (13)  a  constant  loading 
direction  implies  constant  strain  rate  direction  and  thus  constant 

strain  direction.  In  the  combined  tension  a  and  torsion  t  stress  states 

»  ,  •  » 

ojj  is  equal  to  sero  except  On  ■  (2/3)0,  °22  “  °33  *  (”  1/3)0,  and 

°12  “  °21  *  T*  ttultistep  proportional  loadings, 

V  V  V 

t/[(2/3)o]  ■  e12/cll  "  where  k  is  a  constant  for  each  test,  and  Cjj 

and  TniaX  can  be  expressed  as: 

eij  -  {[(3/4)eill2  ♦  (e12)2>1/2  "  ell(k2+(9/16))l/2  »  (17) 

Tmax  “  ((<J/2)2+t2]1/2  -  ((2/3)a] <k2+(9/16))l/2  .  (18) 

v 

By  substituting  (17)  and  (18)  into  (16)  cn  under  a  series  of  m  steps 
in  stress  can  be  obtained  as 

eIl  ■  ( [(2®i/3)  F(Tmax(lpJ  ^  (ti)  +  ... 

1/N 

*  [(20^1/3)  F(Tmax(m-l)]  (t-m-l  ”  ^-2^ 

k/H  N 

♦  [(20m/3)  F(tjmx(o))]  (t-tn-x))  ,  t^x  <  t  .  (19) 

V 

A  similar  form  can  be  obtained  for  e{2  with  20/3  replaced  by  f.  The 
above  expressions  are  exactly  the  same  as  those  reported  in  [10,16] 
with  the  stress  dependent  coefficients  replaced  by  F(  Tmax]°ij* 


DEVELOPMENT  OF  THE  KINEMATIC  HARDENING  FLOW  RULE:  Neglecting  the 
weakening  effects  mentioned  in  a  previous  paragraph,  a  general  form  of 
the  kinematic  hardening  model  can  be  expressed  as: 


®ij  “  [Hl(|lijl)]£ij/|lijl  (flow  rule)  ,  (20) 


and 


®ij  *  H2(l°ijl  ,  leijl)cij  (evolution  rule)  ,  (21) 

•  t 

where  £ij  •  °ij~°ij,  °ij  is  the  stress  deviator,  is  the  state 
variable  (back  stress),  and  laijl,  leijl  ere  the  magnitudes  of 

the  corresponding  stress  tensors,  which  in  principle  can  be  defined  by 
either  Mises  or  Tresca  relations,  see  [6—8] .  By  symmetry  in  tension 
and  compression  and  the  indifference  of  creep  strain  to  the  stress 
direction,  as  implied  by  Eq.  (3)  and  (6),  it  is  reasonable  to  assume 
the  initial  values  of  a^j  are  all  zero. 

I 

In  a  pure  torsion  test,  is  equal  to  sero  except  for 

I  • 

°12  *  °2l  “  T*  which  implies,  according  to  Eq.  (20)  and  (21),  that  all 

•V  »V 

the  tensorial  components  are  sero  except  for  ■  ^21  t  0* 

•  • 

°12  "  *21  *  °»  C12  "  E 21  ^  °»  °12  "  ®2i  *  a  1*  0,  and  that  |ljj|  ,  I j | 
and  I e £ j I  are  equal  to  lT-«l ,  M  ,  and  I *12 1  *  respectively.  Eq.  (20) 
and  (21)  thus  reduce  to  the  following  forms: 

'12  -  Hi(t-a)  ,  (22) 

*  .  V  *V 

a12  "  H2(o,el2)€i2  . 


(23) 
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Since  the  initial  value  of  a  is  zero*  the  initial  shear  creep 

rate  is  a  function  of  shear  stress  only*  according  to  Eq.  (22)*  i.e., 
•V 

e12(t*0)  *  H!(t).  Therefore  function  Hj  can  be  specified  by  the 
relationship  between  the  initial  sheaT  creep  rates  and  shear  stresses. 

However*  the  initial  creep  rates  are  so  rapid  as  to  be 
indetenninant*  possibly  infinite  according  to  Eq.  (1),  (2).  In  order 
to  obtain  a  reasonable  estimation  of  an  initial  creep  rate*  the  first 
two  hours  of  creep  data  of  all  the  tests  reported  in  [13]  were  fitted 
by  Eq.  (1)  with  nj  *  0.223  to  obtain  n2  for  use  in  c^(W  model)  and  by 
Eq.  (2)  to  obtain  N  for  use  in  (VP  model).  In  [13]  Eq.  (2)  was 

o 

used  to  fit  the  first  two  hours  of  creep  data  to  get  c^j  for  the 

determination  of  the  elastic  moduli. 

It  was  found  that  there  existed  small  gaps  between  e£j  and  the 

estimated  elastic  responses  for  some  high  stress  level  creep  tests. 

These  gaps  might  be  due  to  plastic  strain  or  high  strain  rate  effects 

during  loading.  These  gaps  were  neglected  and  no  time  independent 

plastic  strain  was  included  in  the  theoretical  analysis,  see  [13].  To 

compensate  for  these  gaps  by  time  dependent  strain*  a  virtual  data  at 

t  ■  0  calculated  from  E  or  G  was  put  in  each  of  those  creep  data  sets 

whose  effective  stress  magnitudes  were  equal  to  or  greater  than 

137.2  MPa  (19.90  ksi)  as  mentioned  in  the  note  for  Table  3  in  [13]. 

Since  the  deviations  were  small  they  actually  had  very  little  effect 

on  the  data  fittings.  The  averages  of  n2  and  N  are  0.330  and  0.292 

♦V  ♦ 

respectively.  Because  only  the  values  of  and  with  n2  and  N 
fixed  are  involved  in  the  creep  rate  estimations*  the  values  of  n2  and 
N  for  individual  tests  are  not  shown  here. 
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Refitting  the  2-hour  data  by  Eq.  (1)  and  (2)  with  n2  and  N  fixed 
.♦V  .  _  ...  .  ♦ 


yields  the  values  of  ejj  in  Eq.  (1)  and  C£j  in  Eq.  (2)  which  are  shown 
in  Table  1.  Instead  of  the  initial  creep  rate  (infinite)  an  "early" 


value  of  creep  rate  was  used  in  determining  H£.  The  "early"  creep 


rates  were  calculated  either  by  n2eijtj  ^  (or  NC£jti^"^)  or  by 


ctY_"2 


♦  N 


ijti  /t\  (or  eijti/ti),  where  t£  is  a  small  time  such  as  1  sec.  In 

either  case  it  can  be  seen  that  the  ratios  of  creep  rates  for  different 

\ 

tests  depend  only  on  or  eij«  Different  approaches  or  choices  of 
t£  can  only  result  in  a  different  constant  multiplier  for  Hi  whose 
effect  will  be  discussed  in  a  later  section.  The  initial  creep  rates 
based  on  the  first  approach  with  ti  ■  1  sec.  are  shown  in  Table  1. 

The  relationships  between  "initial"  shear  creep  rates  and  T  were 
found  to  be  best  represented  by  the  following  expressions: 

For  the  W  model, 


•V 

el2(t«0) 


C!sinh(T/Ki)  , 


(24) 


or,  for  the  VP  model, 

•V 

c12(t-0)  "  C2sinh(T/K2)  ,  (25) 

where  Ki,  K2  are  19.34  MPa  (2.805  ksi)  and  22.23  MPa  (3.225  ksi) 
respectively;  and  Ci,  C2  are  equal  to  0.06371  and  0.1929  in  percent/h, 
respectively,  see  Fig.  9. 

According  to  Eq.  (20)  and  (22)  the  flow  rules  for  pure  torsion 


loadings  can  then  be  expressed  as 
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e12(t)  -  Cisinh  [(T-aJ/Ki]  ,  (W)  , 


e12(t)  =  C2sinh  [(t-«)/K2)  ,  (VP) 


To  obtain  the  evolution  rule,  an  interrelation  between  a,  a,  and 
•V 

e \2  during  deformation  was  derived  from  constant  stress  shear  creep 
tests.  Here,  e^(vp  Model)  is  used  as  an  example;  the  same  method 
applies  for  eV(W  Model)  except  for  different  constants.  From 
Eq.  (27)  it  follows  that:  for  any  stress  histories, 


•  V  *V  _  1/2 

i-a  ■  K2  ln((ei2/C2)  +  t(Ci2/C2)^  ♦  1}  )  •  (28) 


For  pure  torsion  creep  teats  in  which  T  is  constant,  differentiating 
Eq.  (28)  with  respect  to  time  yields: 


a  *  -K2  e'i2/l(ei2)2  ♦  C?]1^2 


**  V  . 

The  evolution  rule  can  be  obtained  by  expressing  ci2  in  Eq.  (29) 
•V 

in  terms  of  a  and  e^2.  Recalling  that  the  creep  strain  can  be 
expressed  by  e}2  ■  e^2lH  in  the  VP  model,  the  following  results  were 
obtained  by  consecutive  differentiation  with  respect  to  time: 


N  e12  tN_1  , 


N(N-l)  ej2  tN“2 


(31) 


From  (30)  and  (31), 


[(N-l)/t]Ei2 


MX17 1}  /gV  “1/ (N-l)  *V 
(N-l)(ei2N)  Cel2>  e12 


ej2  in  Eq.  (32)  is  equal  to  tF^(t)  according  to  Eq.  (6). 

V  .  *V 

Substituting  t  from  Eq.  (28)  makes  e'12  in  (32)  a  function  of  a  and  t\2 

only.  Introducing  Eq.  (32)  into  Eq.  (29)  yields  the  evolution  rule; 

.  .  v 

i.e.,  the  relation  between  a,  a  and  C12  during  creep.  Even  though  the 
evolution  rule  for  hardening  was  derived  from  creep  tests  at  constant 
stress  it  was  assumed  that  the  same  hardening  nature  would  be 
applicable  for  other  stress  histories. 

For  the  W  model  the  evolution  rule  is  the  same  as  the  VP  model 

+  •  +v 

except  N  is  replaced  by  n2  and  ej2  i*  replaced  by  £{2  ■  TFa( t )  from 
Eq.  (5). 

Since  Eq.  (27)  and  (29)  (with  Eq.  (32)  inserted  in  Eq.  (29))  were 

both  derived  from  constant  stress  creep  tests,  the  kinematic  hardening 

rule,  Eq.  (27)  and  (29),  describe  the  creep  curves  very  well 

regardless  of  values  of  Ci  (or  C2)  employed.  But  different  values  of 

Ci  (or  C2)  result  in  different  values  of  a  during  creep  according  to 

•V 

Eq.  (28).  However,  ^12  *n  E9*  (28)  is  very  small  compared  to  Ci  (or 
C2)  except  at  the  very  beginning  of  a  test.  Thus  differences  in  C}  (or 
C2)  are  generally  unimportant. 


Isotropic  hardening  effects  can  be  taken  into  consideration,  if 
required,  by  developing  an  evolution  rule  for  Ki  (or  K2)  as  in  16,7]. 
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KINEMATIC  HARDENING  RULES  FOR  MULTIAXIAL  STRESS  STATES 

Equation  (27),  (29)  and  (32)  can  be  extended  to  multiaxial  stress 

V 

states  as  follows.  For  eijt  (VP  Model), 


'ij 


c2 C  Zi j / 1 j  t )  sinh(  I  ^  I/K2)  , 


ij 


(33) 


a.  .  = 


(I^I2  *  ^  • 


(54) 


where  bj  *  b2^b(b2)  and 

b2  "  |°i j I  ♦  K2  ^n^lEijl/C2)  +  [< I Gi j  1  /C2)2  +  l]l^2>- 

V 

Similar  expressions  can  be  obtained  for  E£j  (VV  model)  with  N  replaced 
by  n2  and  Ft,  by  Fa,  Eq.  (5).  All  the  magnitude  quantities  are  defined 
by  the  Tresca  relation  in  order  to  be  consistent  with  the  previous 
experimental  observation  that  the  strain  rate  magnitudes  defined  by 
the  Tresca  relation  are  constant  at  any  given  time  for  those  creep 
tests  whose  stresses  are  equal  according  to  the  Tresca  relation. 

The  validity  of  these  multiaxial  relations  can  be  checked  by 
viewing  the  data  and  theory  for  the  first  step  of  each  of  the  pure 
tension  and  the  combined  tension  and  torsion  tests  shown  in  Fig.  1-8. 
The  dashed  and  dash-dot  lines  represent  the  results  predicted  by  the 
above  multiaxial  relations  for  W  and  VP  models,  respectively. 

Since  this  paper  deals  only  with  multistep  proportional  loadings, 

•V  r  .V 

C£j  was  taken  to  be  zero  whenever  Eij  was  smaller  than  zero,  i.e.  c { j 


was  frozen.  More  complicated  stress  histories  will  be  discussed  in 
subsequent  papers. 


PREDICTIONS  BY  THE  THEORIES 

In  summary,  two  theoretical  models  are  presented:  W  and  VP.  In 

each  model  the  time  dependent  nonrecoverable  strain  was  represented  by 

either  the  SH  flow  rule,  Eq.  (13),  or  by  the  KH  flow  rule,  Eq.  (33) 

and  (34).  Eq.  (13)  can  be  integrated  explicitly  for  proportional 

V 

loadings  as  shown  by  Eq.  (19)  for  £11.  The  modified  superposition 

principle  was  used  for  representing  the  viscoelastic  component  under 

VE 

variable  stress  states  as  shown  by  Eq.  (8)  for  A11  IMS1*  subroutine 

DVERK  based  on  the  Runge-Kutta-Veraer  fifth  and  sixth  order  method 
with  automatic  time  step  control  was  used  for  the  numerical 
integration. 

The  results  are  shown  in  Fig.  1-8  in  which  W-SH  are  represented 
by  solid  lines,  VP-SH  by  dotted  lines,  W-KH  by  dashed  lines,  and 
VP-KH  by  dash-dot  lines.  The  back  stresses  “ij  at  the  end  of  each 
step  are  shown  in  Table  2. 


DISCUSSION 

A  comparison  of  the  theoretical  results  with  experimental  data 
shows  that  both  W-SH  and  W-KH  give  very  good  predictions  of  the 
material  responses.  Since  the  nonrecoverable  strains  were  frozen  in 
the  KH  models  whenever  was  smaller  than  zero,  W-SH  is  better  than 
W-KH  in  some  steps  such  as  step  6  of  Fig.  1.  Adding  a  linear 
evolution  rule  for  K|  or  K2  in  the  KH  rules  was  tried,  but  did  not 
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result  in  a  smaller  ajj  at  the  beginning  of  these  steps  as  required  in 
order  for  the  applied  stress  to  induce  further  creep.  This  might  be 
an  inherent  disadvantage  for  the  state  variable  approach. 

Since  the  cVE  component  was  not  included  in  either  VP-SH  or 
VP-KH,  they  could  not  predict  the  recovery  strains  and  the  transient 
creep  strains  upon  reloadings  such  as  observed  in  period  3  of  Fig.  3 
and  4,  etc.  Most  of  the  work  of  other  investigators  based  on  VP 
models  used  a  "static  recovery"  term  (representing  annealing  after 
cold  work)  in  the  evolution  rule  to  simulate  the  transient  strains, 
see  [18] .  The  present  work  shows  that  the  eVE  component  can  account 
for  these  behaviors  without  introducing  annealing  effects  into  the 
analysis. 


CONCLUSIONS 

A  strain  hardening  (SH)  and  a  kinematic  hardening  (KH)  flow  rule 
were  developed  from  constant  stress  creep  and  recovery  data.  Both  are 
in  the  form  of  a  system  of  coupled  ordinary  differential  equations. 
Under  multistep  proportional  loadings,  the  strain  hardening  flow  rules 
can  be  integrated  explicitly  to  obtain  all  the  strain  components  at 
any  given  time.  Associated  with  each  hardening  rule,  data  are 
represented  by  two  theoretical  models:  a  viscous-viscoelastic  (W) 
model  and  a  viscoplastic  (VP)  model.  The  viscoelastic  components  were 
represented  by  the  modified  superposition  principle  for  both  models. 
The  agreements  of  the  W  models  with  the  experimental  observations  are 
generally  very  good.  W-SH  is  better  than  W-KH  in  some  step  down 
stress  states.  The  viscoelastic  component  is  essential  in  accounting 
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for  not  only  the  recovery  strain  but  also  the  transient  creep  strain 
upon  reloadings  and  step  up  loadings. 
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Table  1:  Results  used  for  estimation  of  the  "early"  shear  creep  rates 


TEST 

NO. 

SHEAR 

STRESS , 

♦V 

6 12  * 

(in  Eq.  (1) 

with 

n2  =  0.223, 
n2  =  0.330) 

♦ 

e12* 

(in  Eq.  (2) 

with 

N  *  0.292) 

•V 

Cj2  at  t  *  0, 

(from 

+V  n2-1 
"2*12*1 

t j  =  1  sec) 

at  t  »  0, 

(from 

Nc*  t*1"1 
tj  »  1  sec) 

MPa 

(ksi) 

10-Vh"2 

10_4/hN 

10-2/h 

l0“2/h 

T28 

'*7. 4 

(3.97) 

0.1417 

0.3592 

0.1129 

0.3456 

T29 

+1.4 

(6.00) 

0.3629 

0.6608 

0.2891 

0.6358 

T26 

63.2 

(9.17) 

1.034 

1.699 

0.8237 

1.635 

*??<8 

79.2 

(11.49) 

2.404 

3.514 

1.915 

3.381 

T42 

_ 

101.4 

(14.70) 

7.544 

9.564 

6.009 

9.202 

Table  2:  Applied  stress  deviator  (o'  )  In  each  period  and  back  stress  (a. ,)  at  the  end  of  the  period 


lote:  1.  In  pure  tension  tests,  oij  .  o  except  a22  -  033  -  -(l/2)aU  i  0 

2.  In  pure  torsion  tests,  ay  »  0  except  ai2  “  a21  **  0 

3.  In  combined  teneion  and  torsion  tests,  oij  “  0  except  o22  •  «21  I1  0,  u22  •  “33  ■  -(l/2)*j[  i  0 


FIGURE  CAPTIONS 


Fig.  1:  Test  T23.  X  ■  8  h,  Y  ■  0.022.  Numbers  in  parenthesis  on 

the  ordinate  indicate  the  creep  period  to  which  the  ordinate 
number  applies.  Solid  lines:  Viscous-viscoelastic  with 
strain  hardening;  Dashed  lines:  Viscous-viscoelastic  with 
kinematic  hardening;  Dotted  lines:  Viscoplastic  with  strain 
hardening;  Dot-dashed  lines:  Viscoplastic  with  kinematic 
hardening. 

Fig.  2:  Test  A24.  X  =  8  h,  Y  ■  0.022.  See  caption  of  Fig.  1  for 
key  to  theory  curves. 

Fig.  3:  Test  A25.  X  ■  8  h,  Y  *  0.02Z.  Numbers  in  parenthesis  on 

the  ordinate  indicate  the  creep  period  to  which  the  ordinate 
number  applies.  See  caption  of  Fig.  1  for  key  to  theory 
curves . 

Fig.  4:  Test  T26.  X  -  8  h,  Y  ■  0.02Z.  Numbers  in  parenthesis  on 

the  ordinate  indicate  the  creep  period  to  which  the  ordinate 
number  applies.  See  caption  of  Fig.  1  for  key  to  theory 
curves . 

Fig.  5:  Test  CA31.  X  »  8  h,  Y  ■  0.02Z.  Numbers  in  parenthesis  on 

the  ordinate  indicate  the  creep  period  to  which  the  ordinate 
number  applies.  See  caption  of  Fig.  1  for  key  to  theory 
curves . 

Fig.  6:  Test  CT31.  X  ■  8  h,  Y  -  0.02Z.  Numbers  in  parenthesis  on 

the  ordinate  indicate  the  creep  period  to  which  the  ordinate 
number  applies.  See  caption  of  Fig.  1  for  key  to  theory 
curves . 

Fig.  7:  Tests  A32,  a  tension  creep  test,  and  A37,  a  compression 
creep  test  (designated  by  COMP).  X  ■  8  h,  Y  »  0.02Z. 

Numbers  in  parenthesis  on  the  ordinate  indicate  the  creep 
period  to  which  the  ordinate  number  applies.  See  caption  of 
Fig.  1  for  key  to  theory  curves. 

Fig.  8:  Test  A44,  a  compression  creep  test.  X  ■  8  h,  Y  •  0.04Z. 

Numbers  in  parenthesis  on  the  ordinate  indicate  the  creep 
period  to  which  the  ordinate  number  applies.  See  caption  of 
Fig.  1  for  key  to  theory  curves. 

Fig.  9:  Early  shear  rates  vs.  shear  stresses. 
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